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ABSTRACT 

Observations of the hot gas surrounding Sgr A* and a few other nearby galactic nuclei imply electron and proton 
mean free paths comparable to the gas capture radius: hot accretion likely proceeds under weakly-collisional 
conditions in these systems. As a result, thermal conduction, rather than convection, may be important on 
all scales and affect the global flow properties. The self-similar ADAF solution of Narayan & Yi (1994) is 
generalized to include a saturated form of thermal conduction, as is appropriate for the weakly-collisional 
regime of interest. Conduction provides extra heating and yet it reduces the free-free radiative efficiency of the 
accretion flow (by potentially large factors). These idealized solutions suggest that thermal conduction may be 
an important physical ingredient to understand hot accretion onto dim accreting black holes. Conduction could 
also play a role in reducing the rate at which black holes capture ambient gas and in providing an evaporation 
mechanism for an underlying cold thin disk. 

Subject headings: accretion, accretion disks - conduction - black hole physics - hydrodynamics 



1. INTRODUCTION 

Over the past decade. X-ray observations have made it in- 
creasingly clear that black holes are capable of accreting gas 
under a variety of distinct configurations. There is convinc- 
ing evidence for a hot form of accretion occurring at sub- 
Eddington rates in particular, which contrasts with the clas- 
sical "cold and thin" accretion disk scenario (Shakura & Sun- 
yaev 1973). Hot accretion appears common in the population 
of supermassive black holes in galactic nuclei and during the 
quiescent phases of accretion onto stellar-mass black holes in 
X-ray transients (e.g. Narayan et al. 1998; Lasota et al. 1996; 
Di Matteo et al. 2000; Esin et al. 1997, 2001; Menou et al. 
1999; see Narayan, Mahadevan & Quataert, 1998; Narayan 
2003 for reviews). 

Even today, however, the exact nature, structure and prop- 
erties of these hot accretion flows remain controversial. In- 
spired by the work of Shapiro, Lightman & Eardley (1976), 
Narayan & Yi (1994; 1995a,b) derived self-similai- ADAF so- 
lutions which emphasized the important stabilizing role of ra- 
dial heat advection (see also Ichimaru 1977; Rees et al. 1982; 
Abramowicz et al. 1988). Subsequent analytical work on 
hot accretion flows has emphasized outflows, motivated by a 
positive Bernoulli constant (ADIOS; Blandford & Begelman 
1999), and the potential role of convection (CDAFs; Narayan 
et al. 2000, Quataert & Gruzinov 2000). Hydrodynamical and 
MHD numerical simulations have also greatly contributed to 
the subject by highlighting important dynamical aspects of the 
problem (Hawley et al. 2001; DeVilliers et al. 2003; Igumen- 
shchev et al. 2003). 

The purpose of this Letter is to make the case that ther- 
mal conduction, which has been a largely neglected ingredi- 
ent, could affect the global properties of hot accretion flows 
substantially. In §2, we use existing observational constraints 
on a few nearby galactic nuclei to argue that hot accretion 
is likely to proceed under weakly-collisional conditions in 
these systems, thus implying that thermal conduction could 
be important. In §3, we extend the original 1 -temperature 
self-similar solutions of Narayan & Yi (1994) with a saturated 
form of thermal conduction and study the effects of conduc- 
tion on the flow structure and properties. In §4, we comment 



on some of the limitations of our work, on possibilities for fu- 
ture work and on potential additional consequences of thermal 
conduction for hot accretion. 

2. OBSERVATIONAL CONSTRAINTS AND COLLISIONALITY 

Chandra observations provide tight constraints on the den- 
sity and temperature of gas at or near the Bondi capture radius 
in Sgr A* and several other nearby galactic nuclei. These 
observational constraints have been used before to estimate 
the rate at which the gas is captured by the black hole in 
these systems, following Bondi theory (e.g. Loewenstein et 
al. 2001; Di Matteo et al. 2003; Narayan 2003). Here, we use 
these same constraints (taken from Loewenstein et al. 2001; 
Baganoff et al. 2003; Di Matteo et al. 2003; Ho, Terashima & 
Ulvestad 2003) to calculate mean free paths for the observed 
gas.' Galactic nuclei and the corresponding gas properties on 
1" scales are listed in Tabled ny is the gas number density. 
Till is the temperature, R\ is the 1" size-equivalent at the nu- 
cleus distance, l\ is the mean free path for the observed 1" 
conditions and /^cap is the capture radius, inferred from the 
gas temperature and the black hole mass in each nucleus. We 
have used Spitzer's (1962) expressions for collision times, a 
standard expression for the thermal speed and a Coulomb log- 
arithm In A ~ 20 to recover the simple mean free path scaling 
/ ~ 10'*(rVn) cm of Cowie & McKee (1997; valid for both 
electrons and protons in a 1 -temperature gas). We have used 
this relation to calculate l\ values in Table[2 

The inferred mean free paths are in the few hundredths to 
few tenths of the observed 1" scales {l\ /R\ values in Table^. 
Assuming no change in the gas properties down to the in- 
ferred capture radius, values of /i//?cap > 0.1 are deduced in 
four of the six nuclei. One could also reformulate these scal- 
ings by stating that the 1" mean free paths are systematically 
> 10"*~''/?s, where Rs is the Schwarzschild radius (the typical 
lengthscale associated with the accreting black holes). These 
numbers strongly suggest that accretion will proceed under 
weakly-collisional conditions in these systems. 

' It should be noted that these constraints are not all equally good, in the 
sense that Chandra has probably resolved the gas capture radius in Sgr A* 
and M87 but not in the other nuclei. 
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Table 1. Observational Constraints on Collisionality 



Nucleus 


Tin 




h/Ri h/Rc;p 




(cmr^) (lO' K) 


(cm) 





Sgr A* 


100 


2.3 


1.3 X 10" 


0.4 


0.4 


NGC 1399 


0.3 


0.9 


3.1 X 10^0 


0.009 


0.02 


NGC 4472 


0.2 


0.9 


2.5 X lO^o 


0.016 


0.07 


NGC 4636 


0.07 


0.7 


2.2 X 1020 


0.032 


0.6 


M87 


0.17 


0.9 


2.7 X lO^o 


0.018 


0.02 


M32 


0.07 


0.4 


1.2 X lO" 


0.2 


1.3 



The weakly-collisional nature of ADAFs has been noted 
before (e.g. Mahadevan & Quataert 1997), in the sense of 
coUision times longer than the gas inflow time, but these 
are model-dependent statements. Direct observational con- 
straints on the gas properties near the capture radius favor a 
weakly-collisional regime of hot accretion, more or less in- 
dependently of the exact hot flow structure. Whether the gas 
adopts an ADAF, CDAF or ADIOS type configuration once 
it crosses the capture radius, it is expected to become even 
more weakly-collisional as it approaches the black hole, since 
the relative mean free path scales as l/R oc R~^I^~p in these 
flows with a virial temperature profile and a density profile 
p oc R~^l^^''. It should be noted, however, that tangled mag- 
netic fields provide a way of limiting the weak collisionality 
of the flow, even though their efficiency at doing so is not well 
known. We go back this important issue in §4. 

3. SELF-SIMILAR SOLUTION WITH SATURATED CONDUCTION 

The standard Spitzer formula for thermal conduction ap- 
plies only to gas well into the collisional ("fluid") regime, 
with a mean free path / <C L, for any relevant flow scale 
L. Anticipating temperature variations on local scales ^ R 
from self-similarity, in the weakly-collisional regime of inter- 
est here, with I ^ Ror even R,a saturated (or equiv. "flux- 
limited") thermal conduction formalism must be adopted if 
one is to avoid unphysically large heat fluxes. We adopt the 
formulation of Cowie and McKee (1977) and write the sat- 
urated conduction flux as Fs = 5^spc^^, where 05 is the satu- 
ration constant (presumably < 1), p is the gas mass density 
and Cs is its sound speed. With this prescription, the largest 
achievable flux is approximated as the product of the thermal 
energy density in electrons times their characteristic thermal 
speed (assuming a thermal distribution and equal electron and 
ion temperatures; see Cowie and McKee 1977 for details). 
Because it is a saturated flux, it no longer explicitly depends 
on the magnitude of the temperature gradient but only on the 
direction of this gradient. Heat will flow outward in a hot ac- 
cretion flow with a near-virial temperature profile, hence the 
positive sign adopted for Fs. 

With this simple formulation for conduction, the steady- 
state 1 -temperature self-similar ADAF solution of Narayan 
& Yi (1994) can be generalized to include the diver- 
gence of the saturated conduction flux in the entropy-energy 
equation.^ Adopting the same geometry and notation as 
Narayan & Yi (1994), we consider a "visco-turbulent", 
differentially-rotating hot flow which satisfies the following 

- Honma (1996) and Manmoto et al. (2000) have investigated the role of 
"turbulent" heat transpoil in ADAF-like flows. By relying on a saturated form 
of "microscopic" thermal conduction which is physically well-motivated, our 
analysis differs substantially from theirs. Saturation is key to a self-similar 
scaling, which would not obtain for a standard Spitzer conduction law. 



height-integrated equations for the conservation of momen- 
tum and energy 
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supplemented by the continuity equation, which is rewritten 
as M = -AuRHvp. The last term in Eq is the additional 
divergence of the saturated conduction flux. In these equa- 
tions, R is the cylindrical radius, VIk is the Keplerian angular 
velocity around the central gravity-dominating mass, v is the 
gas radial (in)flow speed, Vl is its angular velocity, p is the 
mass density, c, is the isothermal sound speed, s is the gas 
specific entropy, T is its temperature, H = RcJvk is the flow 
vertical height and / < 1 is the advection parameter (/ = 1 
for negligible gas cooling). A Shakura-Sunyaev prescrip- 
tion has been used to capture the "visco-turbulent" nature of 
the flow, with an equivalent kinematic viscosity coefficient 

= ac^ I fix, where a is the traditional viscosity parameter. 

We look for solutions to these equations of the form 

p = poR-'^\ V = voR-''\ n = noR-'^\ c] = c^R-' , (4) 

and use the notation vk = vkqR'^^^ = R^k = R^koR'^^^- We 
find that the three unknowns vq, flo and Cso satisfy the rela- 
tions 
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while the density scale, po, is fixed by the accretion rate M. 

Contrary to the original ADAF solution, the energy relation 
(Eq. Q) is not a quartic in Csq/vko because of the extra satu- 
rated conduction term. We solve this fourth order polynomial 
with a standard numerical technique (Press et al. 1992) and 
discard the two imaginary roots as well as the real root with 
negative c^o- While in the original ADAF solution the two real 
roots for Cjo were degenerate (in the sense that only mat- 
tered), this degeneracy is lifted up with saturated thermal con- 
duction. The root with negative Cso is physically unacceptable 
because it leads to a conduction flux going up the temperature 
gradient. Note that we have focused on accretion solutions 
(vq < 0) with 1 < 7 < 5/3 in our search, but by analogy with 
Narayan & Yi (1994), one also expect rotating wind solutions 
(vo > 0) with saturated thermal conduction to exist for values 
of 7 > 5/3. 

The main features of the new self-similar solutions with sat- 
urated conduction are illustrated in Fig.^ as a function of the 
saturation constant, We are showing results for two dis- 
tinct solutions here, with 7 = 1.5 and 7 = 1.1, and have fixed 
/ = 1 and a = 0.2 in both cases for definiteness. Standard 
ADAF solutions are recovered at small <i>5 values. As ap- 
proaches unity, however, the solutions deviate substantially 
from the standard ADAF. The sound speed, Cj, and the radial 
inflow speed, -vv, both increase with the magnitude of con- 
duction, while the squared angular velocity, il^, decreases. 
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Fig. 1 . — Changes in the sound speed (Cs/vn, sohd hnes), angular rotation 
velocity (Q^/H^, dashed lines) and radial inflow speed (— iv/i'a:, dotted lines) 
for two distinct solutions, as a function of the saturation constant, <E>,. At low 
<&., values, these solutions match a regular ADAF. The solution extending 
up to log <I> 5 ~ — 0.5 has a gas adiabatic index 7 = 1.1, while the solution 
extending up to log<l?j ~ -1.9 has 7 = 1.5. Both solutions have a = 0.2 
(viscosity parameter) and / = 1 (advection parameter). 



Solutions still exist when Q,- becomes negative but they are 
obviously physically unacceptable. The solution with 7=1.5 
reaches the non-rotating limit at log$j ~ -1.9, while the so- 
lution with 7=1.1 reaches this limit for log'I'j ~ -0.5. Both 
solutions have a sub-virial temperature (Cj < v^) and a sub- 
sonic radial inflow speed (-v^ < Cj), which appear to be rather 
general properties. Supersonic solutions may exist for large a 
values when 7^5/3. 

As the level of saturated conduction is increased, more and 
more heat flows from the hotter, inner regions, resulting in a 
local increase of the gas temperature relative to the original 
ADAF solution. Simultaneously, the gas adjusts its angular 
velocity (which reduces the level of viscous dissipation) and 
increases its inflow speed to conserve its momentum balance 
(increasing at the same time the level of advection). One can 
show from the energy equation that oc (1 - |^cond/?adv|), 
where <7cond and ^adv are the magnitudes of the conduction 
and advection terms: solutions cease to exist once advection 
becomes unable to balance the heating due to saturated con- 
duction. Solutions with = can be seen as Bondi-like with 
saturated conduction. 

We have found that variations in the advection parameter, 
/, have little effect on the solutions, as long as / < 1 . The 

breakdown of the solutions (when — » 0) occurs at lower 
values of the saturation constant $j for smaller a and larger 7 
values. While this breakdown may seem problematic, because 
it occurs for smallish values of when 7 < 5/3 (Fig. 
it may well turn out to be a simple pathological feature of 
the ID height-integrated equations we solved. By analogy 
with the results of Narayan & Yi (1995a) on convection in 2D 
self-similar ADAF solutions, we might expect conduction to 
dominate over advection only in the polar regions of the flow, 
in which case it may be possible to extend the solutions to 




- 



-5 -4 -3 -2 -1 

log[*J 

Fig. 2. — Drop in the level of free-free emission, relative to the standard 
ADAF value, as a function of the saturation constant, for the same two 
solutions as in Fig.fTland a fixed accretion rate, M. 



larger values (and perhaps launch outflows). 

From the solution scalings, it is possible to calculate the 
modified level of emission expected relative to a standard 
ADAF. Assuming that free-free emission is the dominant 
mechanism, it scales as p^T^I^ oc (csoVq)"' for a given M. 
Changes in this quantity as a function of the saturation con- 
stant, $5, for the two same specific solutions of Fig. ^ are 
shown in Fig. |2l (normalized to the low-conduction=ADAF 
values). The emission is effectively reduced, potentially by a 
large factor (specifically by ^ 20 in our 7=1.1 model). Even 
though conduction heats up the gas locally, the reduced den- 
sity resulting from the larger inflow speed dominates, leading 
to a net decrease in the expected level of free-free emission. 
Clearly, by reducing the radiative efficiency at a given M, sat- 
urated conduction could help explain the phenomenology of 
dim accreting black holes. However, work beyond these ide- 
alized self-similar solutions will be required before we can 
elucidate the role of conduction in hot accretion flows. 

4. DISCUSSION AND CONCLUSION 

The possibility that tangled magnetic fields strongly limit 
the efficiency of thermal conduction in a hot flow remains 
a major theoretical uncertainty. In a sense, we have ac- 
counted for this possibility by allowing the value of the sat- 
uration constant, to vary in our solutions. In recent 
years, a similar issue has been discussed in the cooling flow 
context: several studies have concluded that magnetic fields 
would probably not strongly limit thermal conduction (e.g. 
Narayan & Medvedev 2001; Gruzinov 2002; see also Chan- 
dran & Cowley 1998; Chandran et al. 1999), but it remains 
somewhat of an open issue. Explicit numerical simulations 
with anisotropic conduction along fields in an MHD turbulent 
medium could greatly help in settling this issue (see, e.g., Cho 
& Lazarian 2004). 

An important limitation of the self-similar solutions 
with saturated conduction is their 1 -temperature structure. 
Narayan & Yi (1995b) have shown how the 2-temperature 
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property of the hot flow is crucial to obtain self-consistent so- 
lutions with realistic cooling properties. The decoupling of 
the electron and ion temperatures would modify the role of 
saturated conduction in the inner regions of the hot flow, an 
effect which is not captured by our solutions. In this situa- 
tion, the ions may be insulated from electron conduction. In 
addition, self-similarity is broken in the 2-temperature regime 
and the electron temperature profile flattens to a sub-virial 
slope (e.g. Narayan, Mahadevan & Quataert 1998). While 
this suggests a reduced role for conduction, the break-down of 
self-similarity also implies that the heating of the outer flow 
regions must be done at the expense of the hottest inner re- 
gions. Given the strong dependence on temperature of cyclo- 
synchrotron emission, even a moderate drop in temperature 
in the hottest regions could lead to a significant reduction in 
the radio emission properties of the flow. More detailed mod- 
els are needed to address these issues and determine quanti- 
tatively by how much conduction can contributes to the low- 
radiative efficiency of dim accreting black holes. 

In the last few years, much emphasis has been put on reduc- 
ing the rate at which gas in a hot flow is accreted by the black 
hole, either by invoking a positive Bernoulli constant to jus- 
tify the existence of powerful outflows (Blandford & Begel- 
man 1999) or by invoking a net reduction in the accretion rate 
due to convection (Narayan et al. 2000; Quataert & Gruzi- 
nov 2000). The presence of strong conduction in the flow 
may pose a challenge to both of these scenarios: if conduc- 
tion is important on all scales of interest in the hot flow (be- 
cause of the large mean free paths), the Bernoulli constant^ 
becomes irrelevant (as it is a property of adiabatic flows) and 
convection may be suppressed when any displaced fluid ele- 
ment efficiently reaches thermal equilibrium with its environ- 
ment through conduction (but see Balbus 2001). Obviously, 
this is not to say that outflows are not expected from hot accre- 
tion onto black holes, since different configurations than the 
ones described by idealized steady-state analytical solutions 
are certainly possible (as shown explicitly by global turbulent 
MHD simulations). 

If conduction is indeed important for hot accretion onto 
black holes, it could also have other interesting consequences 

Satisfying tiie same Eq |5) as ADAFs (Narayan & Yi 1994), solutions 
with saturated conduction have formally larger positive Bernoulli constant 
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besides the ones illustrated in our specific solutions. For in- 
stance, if the hot flow inside of the capture radius around a 
black hole is able to heat up the ambient gas outside of that 
capture radius, the Bondi rate of gas capture may be effec- 
tively reduced. While Gruzinov (1998) invoked "turbulent" 
heat flows to achieve this, the weakly-collisional conditions 
discussed in §2 for several galactic nuclei suggest that this 
could be done by "microscopic" conduction itself. A reduced 
rate of gas capture below the Bondi value could go a long 
way in explaining the low luminosity of these accreting black 
holes, but here again detailed models are required to address 
this issue quantitatively. Conduction could also have impor- 
tant consequences for the dynamical structure of hot flows, as 
argued recently by Balbus (2001) in a stability analysis in the 
presence of anisotropic conduction. In addition, in a weakly- 
collisional hot flow, it is in principle possible for the "micro- 
scopic" viscous stress tensor to contribute significantly to the 
overall flow dynamics, along with the turbulent component 
usually emphasized in accretion theory (e.g. Quataert et al. 
2002). 

Finally, conduction may offer a solution to the disk evap- 
oration problem. Around both stellar-mass and supermassive 
black holes, there is evidence for inner truncation of thin ac- 
cretion disks (e.g. Esin et al. 2001; Quataert et al. 1999). 
The process which "evaporates" the inner disks into a hot flow 
is not well understood, however. A simple calculation shows 
that, independently of the mass of and distance from the black 
hole, the saturated conduction flux (~ pc,) available to heat 
up a thin disk embedded in a hot flow (using ADAF p and c, 
scalings for simplicity; Narayan et al. 1998) is within ^ a"' 
of the traditional viscous dissipation rate per unit area of the 
disk (for a same M in both components). This suggests that 
saturated conduction is energetically capable of contributing 
to the evaporation of thin disks in both classes of systems. 

The author is grateful to S. Balbus, W. Dorland and 
E. Quataert for useful discussions. This research was sup- 
ported in part by the National Science Foundation under Grant 
No. PHY99-07949 (at KITP). 



because of an increased temperature and reduced angular velocity. 
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